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ABSTRACT 


Supersonic combustion in the inviscid stream has been 
shown to reduce the base drag on a projectile. Fast reacting 
fuel burned externally in the inviscid, supersonic flow gen- 
erates compression waves, which, when properly originated, 
impinge on the recirculation bubble aft of the projectile 
base and increase the base pressure. A free-jet wind tunnel 
was constructed to study external burning where heat addition 
was simulated by nozzle contour. Sufficient increase in base 
pressure was attained to not only reduce base drag but also 
provide thrust. Thus, external burning provides a method 


for increasing the performance of projectiles. 
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INTRODUCTION 


For an axisymmetric object flying at supersonic speed as 
much as 50 percent of the overall drag can be attributed to 
the negative gage pressure acting upon the base. This rela- 
tionship between base drag and the remaining drag on the body, 
which consists of wave, skin friction, and boat-tail drags, 
is observed primarily in the low supersonic flight regime. 
For Mach numbers near 2 the pressure on the base can be re- 
duced to as little as 50 percent that of the ambient air 
while the nose pressure Can be kept relatively low by using 
slender, streamlined shapes, With increasing Mach number, 
however, a limit to the base drag is reached since the base 
pressure cannot be less than a negative one atmosphere, 
whereas the pressure on the nose will continue to increase 
with Mach number, 

The reduction of drag is desirable in order to increase 
the range and terminal velocity of a projectile. Pressure 
or wave drag and skin friction drag have been satisfactorily 
reduced by designing streamlined shapes. Base drag, however, 
has not been reduced significantly although the use of boat- 
tails and base bleed have contributed slightly. The most 
promising method of base drag alleviation is that due to 
external burning. Heat addition in the inviscid flow in 
proximity of the base will create compression waves that, 
if optimally originated, will impinge on the base flow 


region. The interaction of these compression waves will 
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deflect the free streamline toward the symmetric axis, and 
the static pressure in the recirculation zone will increase 


thus reducing the base drag. 
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PEOR 


The base flow problem is visualized most easily with 
the inviscid flow model in Fig. 1 which ignores the effects 
of viscosity. 
A.  INVISCID FLOW 
A semi-infinite, two-dimensional model can be formu- 
lated which easily satisfies all the necessary boundary 
conditions. 
1. ‘Two-Dimensional Model 
If a value of Pr (= base pressure) is selected, 


then the initial turning through the Prandtl Meyer expansion 


is uniquely determined. The pressure and, therefore, the 


EXPANSION FAN 
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BASE 





ΝΕ OF SYMMETRY 


FIG. I 
Two- DIMENSIONAL INVISCID FLOW 
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velocity and Mach number are constant along the streanline 
BC. The Mach number can be calculated from the Prandtl 
Meyer equations. For uniform, two-dimensional flow over a 
convex corner, the pressure depends only on the angle of 
inclination of the free streamline (angle Y in Fig. l), and 
thus BC is a straight line. The shock wave originating at 

C is weakened eventually by the interaction of the expansion 
waves. Near the plane of symmetry, the shock must return the 
flow to the free stream direction. As the flow continues 
downstream of the trailing shock, the static pressure ap- 
proaches the free stream pressure. Unfortunately, for a 
given Μον an infinite number of possible solutions exist. 
As Pp is decreased, the Prandtl Meyer angle is increased, 
and point C moves closer to the base. A limit is reached 
eventually where for a given local Mach number My y the 
trailing shock cannot return the flow to its original 
direction. The maximum deflection shock wave for a given 
local Mach number determines a limiting flow from the infi- 
nite number of possible solutions. This limiting case 
yields the nen possible P, and thus the maximum base 


B 


drag in inviscid flow. 
2. Axisymmetric Model 
Inviscid theory can be extended to axisymmetric 
shapes with some complications. The expansion characteris- 
tics issuing from the base corner are not straight lines as 
in the Prandtl Meyer flow except near the corner. Because 
the flow conditions upstream of the trailing shock are 


dependent not only on the streamline inclination at the 


14 





base corner but also on the history of the flow upstream of 
the Mach lines passing through the corner, the free stream- 
line of constant pressure is not straight, as shown in Fig. 
2. The method of characteristics can be utilized in order 


to determine possible flow patterns. 
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AXISYMMETRIC INVISCID FLOW 


The inviscid flow pattern of an axisymmetric body cannot be 
constructed to the axis of symmetry and still satisfy the 
boundary conditions. This is a consequence of the curvature 
of the free streamline and the fact that for a given Mach 


number, a single trailing shock is limited in the amount of 
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corno Can provide, The angle of inclination of the 
constant pressure free streamline continues to increase 
monotonically as the axis of symmetry is approached. The 
O” pattern, therefore, cannot be continued beyond a point 
where the limiting angle (a in Fig. 2) exceeds the angle 
that the single trailing shock wave can return the flow 
parallel to the axis. In real flows, the wake prevents the 
free streamline from reaching the axis, and thus this limita- 
tion yields some similarity to observed flows. Flow fields 
containing a free streamline which does not reach the axis 
of symmetry can be considered as those that would exist in 
an inviscid flow about an axisymmetric body which has an 
infinitely long cylindrical sting attached at the base as 


Town in Fig. 3. 





FIG. 3 
CHAPMAN ANALOGY BETWEEN WAKE 
DIAMETER AND INVISCID FLOW ROD 
DIAMETER 
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The fraction d/2h equals the ratio of rod diameter to 
body diameter. As d/2h approaches one, the two-dimensional 
result for Pa is obtained. With a rod present, inviscid 
theory predicts a lower base pressure and consequently a 
higher base drag for two-dimensional profiles than for 
axisymmetric shapes. 

This is in qualitative agreement with the experiment. 
Chapman [Ref. 1] found that the base drag of axisymmetric 
bodies in an inviscid flow was so small that it could not 
account for observed drag values determined from experi- 
ment in real fluid flow. This strongly suggested that 


viscous effects were essential in determining base flow. 


B. VISCOUS FLOW 
The inclusion of viscosity effects results in a more 


complicated flow model as shown in Fig. 4. The fluid viscosity 


EXPANSION 


FIG. 4 
VISCOUS FLOW AT THE BASE 
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causes a slow circulatory motion in the dead air region, and 
the viscous mixing process causes the boundary layer to 
thicken as it approaches the trailing shock wave. The base 
pressure arrives at a steady state equilibrium value because 
of the scavenging effect of the outside flow on the recircula- 
tion region. As some of the dead air mass is removed, the 
angle of turning at the base corner is increased, and con- 
sequently the base pressure is decreased. The larger turning 
angle increases the flow velocity outside the boundary layer 
which increases the scavenging action, thereby lowering the 
base pressure. The process is repeated until point C in 
Fig. 4 reaches an equilibrium position. As point C moves 
toward the base, however, the pressure ratio of the trailing 
shock wave increases, making it more difficult for the 
scavenged air and the low velocity air in the boundary layer 
to overcome the pressure rise. An equilibrium position is 
reached which balances these two opposing effects. Viscosity, 
therefore, determines a unique solution from the infinite 
number possible for a given Mach number and Reynoids number 
but, unfortunately, complicates the theory excessively. 
C. CHAPMAN'S CORRELATION 

Viscous theory is too involved to be a practical 
method for estimating base pressure. Instead many correla- 
tions have been attempted to aid in the understanding or to 
create a practical means of predicting base drag. According 
to Chapman's model [Ref. 1] base pressure depends on the 
Mach number and pressure at the base corner, the boundary 


layer thickness yd, the body diameter 2h, and the body angle 
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β at the base. For finite bodies, as shown in Fig. 5, 

profile shape directly affects the Mach number and pressure 
at the base corner since the existence of a bow shock wave, 
for example, would alter the Mach number and pressure prior 


Eowthe corner. 


EXPANSION 


FINETE BODY EFFECTS 





inst 
Mand Pi are needed instead of In and pes Chapman 


offered a simple means of calculating the corrected approach 


Mach number M, and pressure P He suggested extending the 


1 1" 
after body one diameter (2h) and taking the average Mach 
number and pressure along this extension as the corrected 
M4 and Pa. From dimensional analysis, Chapman found that 


CoB =F(M,,6/h, ß ) where CoB is the base pressure coeffi- 


cient and is equal to πη» )}. Viscous effects 
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were included in a Reynolds number dependency of the boundary 
layer thickness ό. Reynolds number, Re = Z Nol JA co , is based 
upon the finite body length L and the free stream conditions. 
les is the free stream air density, Vo is the velocity, and 
Ls the coefficient of viscosity. Assuming a fixed My, 
boundary layer theory predicts boundary layer thickness & is 
proportional to I/NRe for laminar flow or & is proportional 

ES (RIE for turbulent flow. With B= 0, the final result 
from dimensional analysis is Sm = LA ARE?) for laminar flow 
and ES -f(L/h Re ^5 ) for turbulent flow. 

Chapman used data from axisymmetric flow to make a corre- 
lation of M with L/hRe2 for laminar flow, as shown in 
Fig. 6, and L/hRe%) Forgeurbulentstlow, 

The effect of both Re and body shape is apparent, and the 
laminar correlation is reasonably close to the data. For 
turbulent flow, Chapman's correlation again is good; however, 
Re has a much smaller effect on base pressure. This is con- 
sistent with theoretical considerations since |/Re 5 is much 


$ 
less Re dependent than is Re”? ‘ 
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Poe ore to OF SPIN 

The effect of spin, or a body's revolution about its 
axis of symmetry, on base pressure also has been shown to be 
considerable. Unclassified material on this subject is 
lacking, and apparently no purely theoretical approach has 
ever been made successfully. It has been concluded instead 
of analysis, that in order to determine the effects of spin, 
experimental measurements on specific object shapes need to 
be made. 

C. Wieselsberger [Ref. 9] showed in 1927 that the drag 
coefficient of bodies at subsonic speeds could be increased 
by spin. Based on drag data of a projectile and cylinder at 
low subsonic speeds with circumferential velocities up to 
three times the forward velocity, Wieselsberger concluded that 
the increase of total drag resulted from increases in skin 
friction and base drag due to spin induced alterations of the 
boundary layer velocity profile. R. Lehnert [Ref. 9] showed 
that spin influenced the aerodynamic characteristics of pro- 
jectiles at supersonic speeds, but comprehensive information 
concerning the effect of spin was lacking. H. H. Kurzweg 
[Ref. 7) showed a close relationship between boundary layer 
and base pressure, as did Chapman, er prompted a study by 
R. Lehnert and S. Hastings [Ref. 9] to determine the base 
pressure behavior as a functional of spin rate. Experiments 
were conducted at the Naval Ordnance Laboratory in a super- 
sonic aeroballistics tunnel. Cone cylinders of varying 
lengths were spun at rates varying from 0 to 43,000 rpm, and 


base pressure was measured through an orifice located on the 
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sting immediately adjacent to the model base. Spark-schlieren 
photographs were taken in order to study the effect of spin 
rate upon the wake. 

In Fig. 7 base pressure ratio was plotted versus model 


spin rate with the dimensionless model length as a parameter. 
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The upper curve, which corresponds to the shortest model, 

shows a considerable decrease in base pressure ratio with in- 
creasing spin rate. The next lower curve also shows a decrease 
in PR. with increasing spin rate. The three lower curves, 
corresponding to the three longest models, show relatively 
little change in the base pressure ratio with changing spin 
rates. As was determined by both Kurzweg and Chapman, for a 


given Mach number and model configuration, the shape of the 
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curve for base pressure versus Reynolds number is determined 
by the state of the boundary layer. 
Figure 8 shows base pressure ratio plotted against 


Reynolds number for Lehnert and Hastings' data. 
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The Reynolds number was referred to free stream conditions 
and the wetted model length. The first curve shows the zero 


spin data. The large decrease in P, with increasing Re up to 


B 
3.6 million is accompanied by the existence of a laminar 
boundary layer along the entire model length. Further in- 
Creases in Re resulted in a decrease in Por and the transition 
in the wake mixing zone appeared to move upstream. Base 


pressure reached its minimum value when transition occurred 


at the model base and increased only slightly at higher Re 
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when transition occurred upstream of the base. The curves 

for various spin rates show a considerable decrease in base 
pressure with increased spin when laminar flow existed at 

the base. This decrease in base pressure due to spin was 
accompanied by an upstream travel of transition in the mixing 
zone as was the case for zero spin. In the transition region, 
the base pressure is unaffected by spin when the boundary 
layer is turbulent at the base. Increasing spin resulted 

in a slight base pressure increase. The minima of the base 
pressure versus Re curves move toward lower Re with increas- 
ing Spin rate. It can be concluded that spin destabilizes 

the boundary layer and thereby produces premature transition. 
Transition at zero spin was observed to occur at nearly one 
diameter upstream of the base while at 40,000 rpm, tbe transi- 
tion moved to almost two diameters upstream. 

Kurzweg [Ref. 8] and Chapman [Ref. 1] observed that base 
pressure and wake length were correlated.  Decreasing base 
pressure corresponds to decreasing wake length as shown in 
Fig. 9. This was confirmed by Lehnert and Hastings [Ref. 9] 
and suggested a n of measuring base pressure on free 
flight projectiles. Note the similarities of Fig. 9 with 


ENG. 7. 
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E. HEAT ADDITION 

Base drag reduction by external combustion was studied 
Pee. C. Strahle [Ref. 13]. Strahle observed that the injec- 
tion of mass or heat into the viscous recirculation bubble 
resulted in low efficiency since the base pressure rise was 
limited to the free stream atmospheric limit. Combustion in 
the inviscid supersonic stream, however, imposed no such 
limit and, therefore, showed the most promise for reducing 
the base drag or for providing a propulsive mechanism. 

l. Control Volume Analysis 

A control volume analysis [Ref. 6] as shown in 


Fig. 10 gives some intuitive feeling for heat addition in 
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the supersonic stream. Applying the momentum theorem for 
the control volume in the x direction, one has 

CR- RJAe= [PdA - S PdA « f'CP-R)dA «( CP-R)JA (À 
where dA is the projected area perpendicular to the free 
stream direction x. For planar flow and no heat addition, 
pressure P along streamline cd is determined by the local 


flow angle in accordance with the Prandtl Meyer function. 
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Figure 11 shows that the static pressure has a lower value 
in the fast turning case (short recirculation region) as com- 
pared to the slower turning case (long recirculation region). 

me [(P-RISA will be less for a fast turning stream- 


line; and, consequently, Pr also will be lower. Since by the 


continuity equation fPudA = [PudA , only changes in momentum 
S d 
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thickness $** and displacement thickness §* can change the 
momentum flux through de and bc, it is the gage pressure 
integrated over cd and de that largely determines base 


pressure. 
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PRESSURE ON SHEAR LAYER EDGE (REF 6) 


Compression waves emitted from the heat addition zone as 


shown in Fig. 12 increase the static pressure along the edge 
of the shear layer, cd in Fig. 11; and therefore, f" ( P- PR) dA 


is increased thereby increasing P The pressure decreases 


B' 
along cd as the flow turns toward the axis. By adding heat 
in the inviscid region, the shear layer is turned inward 

with nå or less loss of static pressure along cd. Pr is in- 


creased which results in a base drag reduction. 
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2. Strahle Two-Dimensional Model 
Strahle [Ref. 13] restricted his heat addition 
model to two-dimensional, turbulent, supersonic flow, with 
the heat addition zone aligneä with the inviscid strean 
beyond the base corner turn as shown in Fig. 13. 

The heat addition zone is rectangular in shape and is 
assumed far enough removed from the viscous region to pre- 
clude any heated streamlines being entrained in the viscous 
flow region, which would result in a performance loss. 
Strahle chose the mixing theory of Crocco and Lees [Ref. 2] 


memaescribe tne base flow because of the method's success in 


29 





predicting the base pressure. Strahle found that the 
inclusion of external combustion did not alter the base 
flow theory of Crocco and Lees but instead only changed the 
initial conditions for Crocco - Lees governing differential 


equation (Eq. 2.30 in [Ref. 2]). 
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FIG. 13 
STRAHLE MODEL FOR HEAT ADDITION 


With heat addition, the angle of inclination € in Fig. 
14 is no longer a function of we alone Put is also a function 
of both W and X. 

The symbol Wa is the ratio of the x component of velocity 


to the speed of sound a A saddle point singularity exists 


t’ 
through which the trajectory must pass in order that the far 


wake limit can be reached. This forms the uniqueness 
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condition; since without it, an infinite number of solutions 
would exist. This saddle point is called a critical point 

and physically represents the transition point from subsonic 

to supersonic flow in the viscous region. A unique solution 

Br Trcedepyzthe Critical point regardless of the external 
combustion occurring in the inviscid stream. The position of 
the critical point, however, and the influence on the allowable 
initial condition will be altered since & is influenced by 


combustion. 





HEAT ADDITION ZONE 


FIG. 14 X3 ΧΑ 
STRAHLE'S LINEAR THEORY CONFIGURATION X 


Strahle analyzed the combustion region by the linear theory 
of heat addition applied to a supersonic stream. The X and Y 


coordinates defined in Fig. 14 are obtained by rotating the 
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axis through an angle O(o)and are aligned with the rectangular 
heat source. In the uniform field with Mach number equal to 


M the perturbation velocity potential equation is 


B’ 
E. Des t Der = $2 (3-12 /P αρ (2) 


and the pressure coefficient is E 3 -26./Ue, Upon inte- 
grating Eq. (2) where f is non zero only for positive X and 
Y, along the Y = 0 line, the perturbation velocity potential 
is (059 »- VA f f C; Gc xe, de. (3) 


With G(o)unknown and thus the orientation of q , which is the 


H 


heat release per unit volume per unit time, also unknown, the 
Bernoulli equation along a streamline states 

dP/Pur = -du/u = do/ del - dde /u 9 -d(W')/2M [1C -0 Μα] ὦ 
Solving for dø and integrating yields 

9-9, = (Μο - 0(M)+ L27,./Ug (5) 
where Ù is the Prandtl Meyer function. Ó is evaluated at 
O(o) on the Y-0 line regardless of the history of the viscous 
region edge. For large heat addition, this is a source of 
error since the initial direction of the viscous region is 
different from its alignnent at the conclusion-of the τας 
fluence of heat addition. Integrating Eq. (3) and substi- 


Being into Eq. (5) gives 


9-9, = D(My)- D(M) — 9, q0x) (6) 
τε N = (3-DAh¿/B, an AW 
g Cx) = O AS S Ky 

- X-Xi-5, Y; BMC Xx 

= Bab if X, X € Xy 


" BgL Yi b - x «xi ew] RN coe 
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with X FX; t Bg NE 
X, 7 X; t Be(Y;+b) 
X, * Xi * Wt B, Y; 
X4 = X¡ + Weg (Y, + b) 


3. Strahle One-Dimensional Model 
Because the linear theory contains contradictions 
of the assumptions at even modest levels of heat addition, 
Strahle suggested a one-dimensional approach for comparison. 
Figure 15 shows the basic scheme under a small angle assump- 
fron that dA, [dx = ο. -, and ἆ Αο/άκ 5 6Θι- 8 
where Ay, is the area of the heat addition zone and A, is the 


area of the adiabatic free stream. 


FIG. 15 
STRAHLE'S ONE-DIMENSIONAL HEAT ADDITION 
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Assuming the flow on the upper side of Ay "ο. 5) 
contains waves of only one family, then 92361 + Me- OM) : 
Using the one-dimensional flow relations that (dAe/Ae)dx = 
(1 - Me) $ Me (dR /P)dx , Strahle obtains 
9-6, = (ME) - AM) (Ae /PYLU-M M 2]H4R [dx - 4A, /dx. (7) 
dPe/dx is always maintained small by the combustion-affected 
region. Thus, for large heat release, CA ax dominates the 
pressure gradient term. Tf glx) is identified with dA, /dx 
in Eq. (6), then an analogy exists between the one-dimensional 
and linear theory approaches. On the average, dA; /dx = 
b(AV/V)/w where AV/V is the constant pressure volume expan- 
sion of the combustible mixture, b is the heat. addition zone 
capture height, and W is the length of the combustion zone. 
With g(x) assumed in the strongest region where KEXEX, 
in Eq. (6), the equivalence is made that 
AV/V = (3-0) Ah¿/HQ% and Nn.=AV/VWA, (8) 
where Ak, is the heat release per unit mass of fuel, M 
Mene air to fuel ratio, and BAM! Me tclona iy. 
AV/V =f(A, P EE with the most important dependences 
being on M and Tg, which is the base temperature. Strahle 
hoped that by combining linear theory with a one-dimensional 
approach the features of both methods could be included while 
removing most of the objections of each method individually. 

4. Strahle's Performance Parameters 
Strahle computed the force on the base to be 
Fe =2( P- P) h where e 


X = 0) and h is the base half height. The specific impulse 


is the pressure on the base (where 


is defined as Lsp = Fa /me where m; is the fuel flow rate. 
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For small values of force level on the base and letting M-»@ 
Lsp = [24he (3-0/8, Bq bW] (Po/P,) (dP/dt) (9) 
where P = Po/Pr, and |_= a heating parameter, The fact that 
3ο creates an absolute linear heat addition limit, which 
though impractical, does suggest the strong role that the 
ExuEETCcal point plays in the nonlinear regime. 

Figure 16 shows the effect the heat addition zone location 


has on specific impulse. As x; is moved from the base toward 


0.6 







οι 5ο 
Y; = 0.4 





04 FROM REF 13 


0.2 


ο 2 4 6 EE X. 
FIG. I6 

STRAHLES SMALL DISTURBANCE SPECIFIC IMPULSE dP/dn 

VS. COMBUSTION ZONE LOCATION X; 


ENesceritical point, which is located at approximately x = 7.7 
in the linear limit, the specific impulse (dP/dn) reaches a 
maximum when the combustion zone is roughly centered between 
the base and the critical point. A steep drop in effective- 


ness occurs as the zone nears the critical point and extends 
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beyond it. This is expected because any influence upon the 
reattachment process downstream of the critical point where 
the flow is supersonic cannot be transmitted upstream. The 
linear limit predicts that external combustion will result 
in a base pressure rise, the specific impulse is expected to 
be high, and the combustion is wasted if the length of the heat 
addition zone W is too long. 

Strahle recomputed the example in Fig. 16 for finite 
combustion strength, and the resulting axial pressure dis- 


tributions are shown in Fig. 17. 
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FIG. 17 
STRAHLE'S FINITE AXIAL PRESSURE DISTRIBUTIONS 
Four of the curves are for X ΕΤΕΠ corresponds to 


the maximum — location in Fig. 16. The top curve is for 
X: = 0.2. Apparently, for this case, it is best to have X: 


small. The reason for this lies in the nonlinear behavior 
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ΙΓ 1551 pOlnt. in all cases, the critical point be- 
comes entrained in the expansion portion of the combustion- 
affected zone. The combustion length W, therefore, fixes the 
interaction length, and Strahle suspected that the most effi- 
cient placement of the combustion region would be as close to 
the base as possible. He also reasoned that the base pressure 
increase resulted from the entrapment of the subsonic recircula- 
tion bubble by compression waves emitted by the combustion zone. 
The singularity can occur only when the inviscid flow is 
directed inward or for G«QO . The viscous zone thickness, 
hovever, must remain finite. Since the primary effect of 
the combustion region is the downward turning of the flow, the 
last available position for the critical point appears to be 
within the combustion-affected portion of the flow. The 
flow passing the base corner cannot turn as far without com- 
bustion, which in turn indicates a rise in base pressure. 
5. Axisymmetric Heat Addition 
Strahle's two-dimensional heat addition reauires 

modification for axisymmetric application. A method by Fuhs 
[Ref. 6] suggests that the change in base drag resulting from 
external burning is BD, = AP A, where A, is the area of the 
base for both planar and axisymmetric bodies as shown in 
Iu 18. 

Equal amounts of heat are released in both the planar 
and axisymmetric cases. Heat addition in the planar case is 
represented by the lines above and below the base area and 
in the axisymmetric case by a concentric circle of radius ro. 


Since in both cases, the length of heat addition zone is 


2) 





zn, cic pase area influenced by heat release is in the 


ratio 


axisymmetric base area rr? . 
planar base area i A n = n / 2. N, (10) 


ΕΝ HEAT ADDITION ZONES 
[| | 


to 
7 


FIG. 18 


COMPARISON OF BASE AREA FOR PLANAR AND 
AXISYMMETRIC GEOMETRY 


The axisymmetric geometry apparently has less base area 
influenced by heat addition. Flow through an annular heat 
addition region was investigated. It was found that the 
compression waves merged sooner with axial symmetry, and that 
in contrast with the planar case, there is an optimal r/r) 


το 1ο. 
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III. APPARATUS 


In order to study the effects of heat addition on the base 
pressure of an axisymmetric body, a free jet wind tunnel was 
constructed at the Naval Postgraduate School's blow down 
laboratory. The test model simulated a 5-inch projectile at 
23,000 feet and a Mach number of 2.0, which are typical condi- 
tions for a maximum range trajectory. Similarity parameters 
were matched by using a Mach 2.0 free jet coaxial nozzle 
shown in Fig. 19 which gave a Reynolds number based on radius 


EN x10. 


RICAJ 
COAXIAL CYLINDER/ NOZZLE 





Do 





A turbulent boundary layer was assumed, and the displace- 
ment thickness 4% was computed to be 0.04] inches at the 
projectile base. This provides a similarity parameter of 
displacement thickness to radius §*/h of 0.0164. 

The equipment consists basically of a 300 psig air supply, 
a glcbe valve, a pressure reducer valve, a 90° elbow, a center 
cylinder support/flow straightening section, a nozzle and an 
overhead exhaust hood. A 4-inch pipe feeds the air from the 
supply tank to the on/off globe valve. The pressure reducer 
valve is followed by a 90? elbow which directs the air flow 
vertically and provides a feedback pressure for the regula- 
tion of the pressure reducer valve. A 2-inch spacer is 
located between the elbow and the center cylinder support 
section and provides internal access to the centerline cylinder. 
The centerline cylinder is supported inside the center cylinder 
support section by flow straighteners located at the top and 
bottom of the 3-foot section of 4-inch pipe. A section of 6- 
inch pipe houses the nozzle on top of the cylinder support 
section and provides armor protection and a support for the 
pressure probe as shown in Fig. 20. Figure 21 shows the 
exhaust hood above the apparatus, the mercury manometer board 
used to obtain pressure data, the mercury vapor light source 
and one mirror used in the schlieren system for flow 
visualization. 

Figure 22 shows the cylinder extending from the cylinder 
support section with the nozzle and armor section removed, 
along with the pressure gage connected to the air supply 


reservoir. 
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TEST APPARATUS 
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The cylinder base is formed by a cap which contains six 
pressure taps. Four taps are located across the base radius, 
and two taps are located just upstream of the base corner. 

The pressure from each tap is carried downward through 
the cylinder and the cylinder support in plastic tubes which 
exit the apparatus through the 2-inch spacer located on top of 
the elbow. The tubes are connected to a mercury manometer from 
which the pressure data was read. Total and static pressure 
also can be measured by the external pressure probe, which 
will be the primary base pressure source when cylinder spin 
is included. The results of this study were for the static 
case with all similarity parameters matched except for spin. 
However, because spin appears to have a considerable impor- 
tance on base pressure, provisions were included to allow 
Spinning the projectile. When a 5-inch projectile is fired, 
the barrel rifling spins the projectile to 250 revolutions 
per second. This establishes a spin similarity parameter of 
ræðu = 0.1595 where r is the radius of the projectile, Wis 
the angular velocity and wis the free stream velocity at 
Mach 2.0 at 23,000 feet altitude. The similarity parameter 
will be matched in the free jet experiment by spinning the 
cylinder at 48,000 rpm. The spinning cylinder design con- 
sists of a 6.5-inch stainless steel cylinder which slides 
over a support cylinder containing high pressure nitrogen 
orifices to provide gas bearing support and align the 
Spinning cylinder radially. A small internal impulse turbine 


driven by high pressure nitrogen will spin the cylinder. 
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Heat addition was simulated by nozzle contour. As de- 
scribed in Ref. 6, the external burning process can be 
simulated by nozzle wall contours vhich generate character- 
istics identical to the supersonic combustion of an external 
burning assisted projectile. A gas generator within the 
projectile injects fuel-rich, hot, combustible mixtures into 
the supersonic inviscid flow where ignition occurs due to 
high temperature and the chemical nature of the fuel. The 
region of heat addition was assumed axisymmetric about the 
cylindrical body of the projectile with the radius τ, in 
Fig. 18 equal to the body diameter. For stoichiometric fuel 
mixtures, calculations showed that the combustion would result 
in a 4.5? deflection of the streamline bounding the heat 
addition zone, and it was this 4.5? deflection inward of the 
streamline at r, from the cylinder centerline which governed 
the wall contour. Half stoichiometric nozzles also are 
planned for future analysis. An iterative procedure using 
the method of characteristics was used to design the nozzles 
to give the proper -4.5° flow angle. Using W to denote the 
length of the combustion zone (see Fig. 14) and h, the cylinder 
radius, several nozzles of varying W/h were fabricated with an 
epoxy cast ES over a wax plug to the correct contour. 
These nozzles provide several data points shown in Table I 
which show the effect external heating has on base pressure. 
During testing, spacers were used to vary the axial location 
of the heat addition zone with respect to the e corner. 
Figure 23 is a photograph of the nozzles, spacers, and tunnel 


armor section. 
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The symbol l represents the distance from the corner 
where the first fully developed compression wave impinges 
the cylinder wall or the wake downstream of the base, and h 
again is the cylinder radius as shown at the top of Figs. 
40, 41, and 42. Five 1/h values were observed for each 
nozzle as shown in Table I. Negative 1/h denotes compression 
beginning upstream of the base corner, positive 1/h signifies 
compression beginning downstream of the corner, and 0 denotes 
compression starting at the corner. Table I is the test 
matrix obtained to date showing the results of the completed 


tests. 
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IV. IEST PROCEDURE 


After the proper nozzle and spacer combination was 
installed on the apparatus and the air reservoir was charged 
to 300 psig, a data point for the test matrix was taken. The 
procedure involved opening the globe valve and setting the 
proper output pressure from the pressure reducer valve. A 
pressure tap in the elbow was connected to a separate 
pressure gage which was used to set the conditions for each 
run. This tap also provided feedback pressure for the 
pressure reducer valve. With the proper output pressure set, 
a schlieren concave mirror system was used to photograph 
each test run as shown in Figs. 25-39, Polaroid positive/ 
negative film was used in order to facilitate print 
reproduction. Base pressure information was read from a 
mercury manometer and recorded for data reduction. At the 
conclusion of the data readings, the pressure reducer valve 
was turned off, and the air reservoir was recharged for the 


DOC run. 
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V. RESULTS AND CONCLUSIONS 


The pressure ratios (Pr/P LG) as shown in Table I indicate 
the performance of the observed simulated external burning. 


Becure 24 is a plot of the results, 
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The data points taken for the two uniform flow nozzles 
compare favorably with M = 2 results compiled by Przirembel 
and Page [Ref. 10]. These points represent the base pressure 
ratio of the cylinder with no heat addition; and because 
Pp Poo is less than 1.0, base drag results. All other data 
points, except one, show an increase in base pressure ratio 


with the applied combustion, indicating a reduction in base 
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Greg. the data for Wh = 4.34, and the last two points 
for W/h = 1.84 in fact, indicate thrust developed at the base 
since P5/ Pg for these cases is greater than 1.0. The data 
indicate, as expected, that a longer heat addition zone 
(Larger W) will result in a higher PIP increase, since a 
longer interaction length of compression waves acting upon 

the recirculation bubble has more effect than a shorter zone. 
It also appears that maximum benefit from the heat addition 
occurs when the entire recirculation bubble is embedded in 

the combustion-created compression waves. For W/h = 4,34, 

a decline in performance is noted as 1/h exceeds -1.0, which 
indicates that for a long heat addition zone, optimal axial 
placement is necessary so as to not waste compression waves 
downstream of the sonic point. Because supersonic flow exists 
downstream of the sonic point, any compression effect occurring 
beyond that point will not be transmitted upstream into the 
subsonic circulation bubble, which would be necessary for a 
base pressure increase. For the shorter heat addition zones 
(W/h = 1.84 and 0.72), the base pressure ratio increases as 
the zone is moved rearward. It is assumed that no compression 
was wasted beyond the sonic point up to 1/h = 0.5. Although 
not tested, a decline in performance is expected with further 
rearward movement of the heat addition region with respect to 
the cylinder base corner. Figures 25-39 show a schlieren- 
derived photograph of each data point obtained with a 
Pretorial description of the heat distribution actually 
Simulated. The double cross-hatched regions depict anomalous 
compression, and the single cross-hatched areas show the heat 


addition zones without multiple reflection. 
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FIG. 25 
PHOTOGRAPH AND SCHEMATIC FOR W/h= 4.34, J/hz -223 
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FIG.26 
PHOTOGRAPH AND SCHEMATIC FOR W/h= 4.34, 1/h=-0.95 
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FIG. 27 
PHOTOGRAPH AND SCHEMATIC FOR W/h= 4.34, 2/h= -0.50 
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FIG. 28 
PHOTOGRAPH AND SCHEMATIC FOR W/h=4.34, 4/h=0.0 
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FIG. 29 
PHOTOGRAPH AND SCHEMATIC FOR W/h= 4.34, Y/h= 0.50 
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FIG. 30 
PHOTOGRAPH AND SCHEMATIC OF W/h=1.84, X/h=-2.23 
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FIG. 31 
PHOTOGRAPH AND SCHEMATIC FOR w/h=1.84, 2/h=-0.95 
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FIG. 32 
PHOTOGRAPH AND SCHEMATIC FOR W/h=1.84, 2/k:-0.50 
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FIG. 33 
PHOTOGRAPH ANO SCHEMATIC FOR W/h=1.84, 2/h=0.0 
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FIG. 35 
PHOTOGRAPH AND SCHEMATIC FOR W/h=0.72, 2/h=-2.23 
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FIG. 36 
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FIG. 37 
PHOTOGRAPH AND SCHEMATIC FOR W/h=0.72,2/h=-0.50 
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FIG. 38 
PHOTOGRAPH AND SCHEMATIC FOR W/h=0.72, #/h= 90.0 
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FIG. 39 
PHOTOGRAPH AND SCHEMATIC FOR W/h=0.72, 2/h=0.50 
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A study of the pictures shows that a uniformly distrib- 
uted heat source was not obtained as desired. Instead, the 
starting compression waves, which originated immediately 
downstream of the nozzle throat, reflected from the cylinder 
wall to the nozzle corner and then into the base flow region 
as shownat the top of Figs. 40, 41, and 42. For the W/h - 
4.34 nozzle, the simulated heating distribution is shown in 
the middle of Fig. 40. 

Compression, originating from the simulated heat addition 
zone, starts at point b and increases to uniform strength at 
point e where it levels at a constant value until just prior 
to point h. The reflected compression waves from the nozzle 
corner superimposed upon the uniform heat cause a spike which 
is approximately twice the magnitude of the intended uniform 
value. This spike is seen in the W/h = 4.34 photographs as 
a shock wave, formed by a coalescing of the reflected 
compression waves, as seen in Figs. 25-29. 

Figures 41 and 42 show the expected compression wave re- 
BIGetion pattern and heat distribution for W/h = 1.84 and 
W/h = 0.72 respectively. Basically, the same phenomena of 
compression wave reflection occurred. The shorter heat addi- 
tion zones, however, ended prior to the reflection, and, 
therefore, the reflected waves created an almost equal 
strength combustion region downstream of the intended heat 
addition zone rather than a spike as seen for the W/h = 4.34 


case. 
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FIG, 40 
COMPRESSION WAVE BEHAVIOR, HEAT DISTRIBUTION AND 
CRITICAL POINT LOCATION FOR W/h=4.34 NOZZLE 
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FIG. 42 
COMPRESSION WAVE BEHAVIOR, HEAT DISTRIBUTION AND 
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The photographs show that for the W/h = 1.84 case, the 
reflected compression waves never coalesce into a shock but 
are spread over a slightly longer axial distance. For W/h - 
0.72, a shock did form, but its strength is notably less than 
the shock occurring for W/h = 4,34. 

Although the reflected compression waves existed for each 
data point, it is assumed that the base pressure was affected 
by the reflection for only the points circled in Fig. 24. 

The circled points were determined by the location where the 
reflected wave impinged the recirculation bubble. As pre- 
viously described, base pressure can be influenced by 
compression waves only if the waves interact with the base 
flow prior to the critical point. Any interaction downstream 
of the critical point will not influence pressure at the base. 
The bottom of Figs. 40, 41, and 42 shows the estimated 
critical point location for the three W/h nozzles. 

A vertical line extended through the intersection of a 
given 1/h value and the critical point curve will have on 
its left the portion of the heat addition region affecting 
the base pressure. 

In conclusion, it is evident that external burning is 
capable of not only reducing base drag but also producing 
thrust under the proper conditions. The optimal benefit is 
derived when the entire recirculation bubble is embedded in 
combustion created compression waves. Waves beginning up- 
stream of the base corner have limited value for raising base 
pressure, and waves interacting downstream of the critical 


point are wasted. 
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